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Theory for a gridded inertial electrostatic confinement (IEC) fusion system is presented that 
shows a net energy gain is possible if the grid is magnetically shielded from ion impact. A simplified 
grid geometry is studied, consisting of two negatively-biased coaxial current-carrying rings, oriented 
such that their opposing magnetic fields produce a spindle cusp. Our analysis indicates that better 
than break-even performance is possible even in a deuterium-deuterium system at bench-top scales. 

The proposed device has the unusual property that it can avoid both the cusp losses of traditional 
magnetic fusion systems and the grid losses of traditional IEC configurations. 


I. INTRODUCTION 

Inertial Electrostatic Confinement (IEC) is a method 
of producing nuclear fusion by trapping and heating 
ions in a spherical or quasi-spherical electrostatic po- 
tential well m- Various schemes have been imple- 
mented 0 0 0 ED3 [IT] , the most common of which - 
known as gridded IEC - employs a concentric spherical 
cathode and anode arrangement as shown in Fig. [l] The 
inner electrode, shown in blue, is the cathode and must 
be highly transparent to converging ions. 

Although one can apply a sufficiently large electric field 
to accelerate ions to fusion-relevant energies, many recir- 
culating ions are intercepted by the cathode grid, result- 
ing in significant heating loss. At bench-top scale, the 
losses are at least five orders of magnitude larger than 
the fusion power produced, even when the devices are op- 
erated in the so-called “star-mode”. Nevertheless, these 
devices have found application as neutron sources HUD- 



FIG . 1: A schematic diagram of a spherical gridded IEC de- 
vice, where the inner electrode is the cathode. 
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Other schemes have been proposed nans] in attempts 
to circumvent grid losses. Such methods aim to replace 
the grid with a virtual cathode formed through the con- 
finement of electrons. Specifically, applying a positive 
bias to the inner grid results in a central electrons fo- 
cus and consequent virtual cathode. A potential well 
is subsequently formed between the anode grid and the 
electron-rich core. Although relatively deep potential 
wells have been formed and significant ion heating ob- 
tained in this way, electron loss to the grid limits the 
attainable efficiency. 

Alternatively, a virtual cathode can be formed through 
the trapping of electrons in a three dimensional magnetic 
cusp system, known as a Polywell H3HSI. This consists 
of an arrangement of current loops that are positioned 
on the faces of a regular polyhedron, most commonly 
a cube. The currents are applied so that a magnetic 
null is obtained at the center of the device. Electrons 
introduced into this device are trapped by the magnetic 
cusp configuration for sufficient time as to form a virtual 
cathode. Inefficiencies arise when electrons are scattered 
into the loss cones of the face, edge, and corner cusps, 
allowing them to escape the trap. 

In this paper, we present a new approach to the reali- 
sation of nuclear fusion in IEC devices. We combine IEC 
and magnetic cusp confinement to yield a hybrid electro- 
magnetic reactor. Ion confinement and heating is pro- 
vided through the electric field, whilst the magnetic field 
serves both to enhance electron confinement and to shield 
the (cathode) grid from ion impact. Earlier attempts to 
combine magnetic and electrostatic confinement [37] did 
not operate in this fashion. 

In the following sections we derive expressions for en- 
ergy gain and loss mechanisms for such a magnetically- 
shielded-grid (MSG) fusion system. Energy gain comes 
from fusion, whilst particle losses through magnetic 
cusps, interspecies energy transfer, ion and electron 
losses, and bremsstrahlung contribute to energy loss. 

Notwithstanding a discouraging prediction of ex- 
tremely poor performance from an early exploratory pa- 
per by Rider [33] , significant net energy gain is predicted 
for a judicious choice of the potential applied to the grid, 
the current in the field coils and the energy of the inter- 
acting species. This is true even for small-scale devices 
and there is little reason to presume it would not also hold 
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true for exotic fuels such as proton-Boron or Helium-3, 
as the Z 2 increase in Bremsstrahlung for these reactions 
would be compensated to some degree by the increased 
energy per fusion event. 


II. DERIVATION OF ENERGY GAIN AND 
LOSS EXPRESSIONS 

A. Device geometry and operation 


practical considerations. 

In general, the device is similar to a spindle cusp ar- 
rangement but with a negative voltage bias imposed on 
the rings. Although other grid configurations can be 
used, the arrangement shown here has cylindrical sym- 
metry that makes analysis and construction relatively 
straightforward. 

In the following subsections we will quantify the var- 
ious energy gain and loss mechanisms, then finally inte- 
grate all of them to illustrate a net fusion energy gain. 


We will examine the simple grid arrangement con- 
sisting of two coaxial negatively biased rings within a 
grounded cylindrical vacuum chamber as shown in Figj2j 
which has been shown PE] to give IEC plasma prop- 
erties. This grid replaces the cathode shown in Fig(l] 



FIG. 2: Two-ring cathode grid with embedded magnetic field 
coils that produce a magnetic cusp. Magnetic field lines are in 
black; electric potential is shown as a heatmap with blue rep- 
resenting a negative potential and red representing a ground 
(a grounded cylindrical chamber is implicit at the edges of the 
diagram) 


Current-carrying coils are embedded within each ring 
and electrically isolated from it. Current circulates in op- 
posite directions in each coil so as to produce opposing 
magnetic fields. The magnetic field created by the coils 
will be shown to shield the surface of the rings from ion 
impact. Moreover, the opposing currents create a mag- 
netic null at the midpoint of the axis through the rings, 
point cusps at the center of each ring, and a line cusp in 
a plane that bisects the axis at right angles. 

The purpose of this magnetic configuration is to in- 
crease the ion density beyond the limit imposed by a 
concentration of only positive charges. Trapping elec- 
trons in this cusped geometry, such that their density is 
lower than those of the ions, enables the ion density to 
be increased while maintaining the net positive charge 
on the non-neutral plasma. This increase in ion density 
will enable an increase in fusion density. Although elec- 
trons are not needed in order to show a net energy gain, 
they enable an increase in the net energy produced for 


B. Density profiles and fusion power 

The conserved quantities (constants of the motion) for 
our system are the total energy, H , and the z— com- 
ponent of angular momentum, L z , which correspond to 
time- and rotation-invariance, respectively. In cylindrical 
polar coordinates (r, #, z) with z the axis of symmetry, we 
have: 


H = t m(v l + v 2 z + v 2 e ) + q</> (1) 

L z = r(mvg + qAg), (2) 

where (j) is the electric potential and Aq is the only non- 
vanishing component of the magnetic vector potential 
for our system. Density variation emerges through the 
position-dependence of Aq and <fi. L z is seen to suppress 
density around the r = 0 axis. 

Where the magnetic field is slowly- changing over dis- 
tance scales of order of the gyroradius, we have a third 
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(“adiabatic”) invariant: /i ~ v ± is the component 

of particle velocity perpendicular to the local magnetic 
field. Adiabaticity will tend to reduce cusp losses (38]; we 
will assume for the sake of conservatism that our plasma 
is nonadiabatic. 

Since any function whose only arguments are constants 
of the motion is a stationary solution to the Vlasov equa- 
tion, we can write electron and ion distribution functions 
in the form / = f(H,L z ) if we neglect individual col- 
lisons and treat only collective behaviour. The collisional 
effects will later be estimated at the steady state of the 
collisonless system. 

For this paper we consider pseudo-Maxwellian distri- 
butions of the form 


f(H,L z ) =n 0 


2m 

2t r 3 fcT 


exp{— H/kT} exp {—L 2 /(2mrQkT)} 


x erfc (/3(H — H c )/kT) 

(3) 

Here ro represents a characteristic radius at which an- 
gular momentum becomes significant; one might expect 
this to be of the order of the particle gyroradius. Spatial 
dependence is also implied through and Aq and thus H. 
This form for the distribution, though very simple, will 
suffice to demonstrate the parameter-dependence of our 
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system. The gradient and value of the cutoff are deter- 
mined by /? and H c respectively. Beyond this cutoff, par- 
ticles are assumed to be able to escape from confinement. 
These pseudo-Maxwellian distributions can be expected 
to yield an underestimate for achievable fusion rate. 

Recognising the symmetry between v r and v z in Eq. [3j 
we can rewrite it as 


2 n 0 / 

m 

y/n V 2kT VL 


l 2kT • 

x exp < 

x exp <j 

r 1 , 

[ 2 mkT 

x erfc I 

( /3m 2 



2 i „.2 


2 qcp 

m 


( r/r 0 ) 2 (mvg + qAgf 


(4) 


with v\ = v 2 + v 2 . In this ( vl,vq ) velocity-space, the 
angle-averaged velocity difference between two points Vi 
and v 2 is 

(vi - v 2 ) = \Jv\ L +vl L v L + (vie ~ v 2 g)vg (5) 
Changing scales to remove dimensions, let: 


v v' = v/y/2 kT/m 

Aq — )> A' e = eAo / y/2 mkT 
r r' = r /tq 


Dropping the primes (we will now work in the rescaled 
quantities): 


at the approximate form, derived in appendix A: 
limn(A g ,(j),r) = ^ exp j- (rZA g ) 2 } j 

\ exp {(r 4 /r 2 )Z 2 Aj — Z</>} C( T -rZAo 1 rw) 

exp {r 2 Z 2 A 2 e — H c } C(rZAo,rw) j> 

where w = 5ft {^J H c — Zcj^j >0, r = r 2 + 1 
and C(a, b ) = erf (a + b) — erf (a — b) 

(8) 

Imaginary values for w are forbidden as they would 
represent plasma with potential energy beyond the ener- 
getic cutoff at H c . 

A sufficient condition for magnetic shielding is C 1 
near the coil surface, i.e: 


Z Aq > \ W - 1 

r 


(9) 


. Eq. [9] is more readily satisfied as the device becomes 
larger. 

Restoring dimensions, we can evaluate Eq. [9] along the 
plane of a coil at half the coil radius, using the approxi- 
mation for a current loop of radius a [39] : 


A 0 (r = a/2, z = 0) « (A^) Mo I 


( 10 ) 


After a small amount of algebra, Eq. [9] reduces to the 
following (in S.I.): 


I> 


29 v / 2m (yH c - (/</ co n + 2VkT (^)) 
4mo<? 


( 11 ) 


f v (VL,V 0 ,Ag,(j),r) 


z 'n) r / 2 I 2 

—j=v L exp \ — {v L + v e 
v 7r 


+ Zd)} 


x exp {-r 2 {v e + ZAg) 2 } 
x erfc (/3 {v 2 L +vj + Z<j>- H c }) 

(6) 

with Z the atomic number (1 for a deuteron, —1 for 
an electron). 

In the limit where H c oo, the density is given by 


n(A e ,(j),r) 



lim f v (v L , vg, Ag, (f>,r) 
H c — yoo 


n 0 


\Jr 2 + 1 


exp 


{- 


+ z<$> 


■)} 


(7) 

The presence of the cutoff complicates this somewhat; 
after some algebra and taking the limit /3 1, we arrive 


For a system of two opposed coils spaced approxi- 
mately a radius apart, the current required differs from 
this by only a few percent. It is worth noting the scaling 
with the square root of energy. 

Given the ion distribution /*, fusion power is given by 
the standard expression 

Pius = I-Efus J d 3 x rii(x) 2 (av) , (12) 

where Ef us is the fusion energy and 
(cry) = (1 /rii(x) 2 )x 

II dv II c^v, IK v_v, )ll CT (IK v_v, )ll)^( v )^( v, )• 

(13) 

To qauntify the above we adopt the fit to the D-D 
fusion cross-section due to Duane [36] . 
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FIG. 3: Example Vacuum electric potential of charged cusp. 


between fields F> 0 , 0o and B' , 0 7 . Let 

Vo = V_L 0 +V || 0 (15) 

v 7 = v 7 ± (16) 

where v 0 and v 7 , respectively, are the velocities at Bq 
and B ' . 

Energy and magnetic moment conservation, respec- 
tively, yield 


\ mv '_ l 2 + Q<f>' = \mv ^ + q<\> 0 

(17) 

2 / 2 


'GLO v' ± 

Bo B' 

(18) 

Combining and rearranging, we get 


Bo vi_l 

B'~ vi + ^o-n) 

(19) 



X < cm > Z (cm) 


which yields the corrected loss-cone angle as 

sin 6 rn = — — (l + — -y ( 0 o — 0 m)^ = - 5 — ( 20 ) 

-Dm \ IILUq j n rn 

where 6 is the angle between the velocity vector v 0 and 
the magnetic field in the low-field region, and L> m is the 
maximum magnetic field strength encountered by the 
particle during its motion. If we can hold the coils at 
a negative potential with respect to the centre of the de- 
vice, the combined electromagnetic trapping is sufficient 
to eliminate electron cusp losses for electrons of energy 
below a small multiple of the virtual anode height. De- 
pending on the geometry, such a virtual anode can form 
naturally when the coils are held at a negative potential 
with respect to the boundary (e.g vacuum chamber). 


FIG. 4: Electric potential of example system with plasma 

present. 


C. Cusp Losses 

Since typical operation of an electromagnetic cusp de- 
vice involves the formation of a central potential, we need 
to obtain the appropriate expressions for cusp confine- 
ment where an electric field is also present. For clarity, 
we work in S.I. in this section. 

Derivation of the loss-cone angle begins with the force 
on a diamagnetic particle: 


D. Grid Losses 


Since the magnetic field of a cusp is convex inward, 
both the curvature and \7B drifts are nowhere towards 
the coils. We can prevent a particle from reaching the 
coils by insisting that the gyroradius becomes smaller 
than the distance of approach. 

We recall the definition of the gyroradius 


r 
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mv 

~qB 


( 21 ) 


Approximating the magnetic field near the current- 
loops with that due to a line-segment of current, we have 


F || = -/iV5 - (14) 

Since charge is conserved, energy conservation yields 
the separate conservation of the magnetic moment fi. 
We consider the motion of a particle turning around 


B(s) 


im9 

2l TS 


( 22 ) 


with s the perpendicular distance from the line- 
segment. 
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We require r g (r) < r — r c for some r c representing the 
thickness (minor radius) of the coils; solving for I yields: 

^ 27TTOV_L _ 2'K^/2E^m i 

Mo9(l - m(i - r v) 

Inspecting Eq. [23j we see that the required number 
of Ampere-turns depends on the achievable separation 
between metal surfaces. This suggests the use of coils 
with a small minor radius. 

If we take r' = 2 r c with r' the distance at which our 
condition applies, then the above simplifies to 


j 4 TTy/2E^mj 

Mo Qi 


(24) 


We can rearrange Eq. 24 to find a critical energy for 
a given current, below which particles cannot escape to 
the grid: 


= Mpgp 2 

327 T 2 rrii 


(25) 


The structure of this agrees with our earlier approach 
through the distribution function (Eq. ED- 

If E cr i t > —qi<p coin upscattered ions will be lost to the 
chamber walls before they are lost to the grid. 


E. Bremsstrahlung 


F. Collisional up-scattering 


Interspecies energy transfer and collisionally-induced 
upscattering both represent energy-loss mechanisms for 
an IEC-style device. We can obtain useful estimates of 
the magnitudes of such losses from the Landau equation 
for the scattering of particles of distribution f a due to 
distribution f b [35] : 


^ = -Vv-J(/ 0 ,A,v) 

at m a 

with 

J (/a,/b,v) = 


(l_ 

uu\ 

V u 

u 3 J 


/b( v ')v 7 f A( V )v 7 

V v/a v v 


m a 


m b 


(28) 

where X denotes the unit dyadic, uu is the tensor prod- 
uct of u with itself, and 


Cab — 


ZlZle^ ln(A) 


Sir c 


2 

0 


(29) 


ln(A) = ln(47re 0 A:T— (v a 2 )\ D /q 2 e ) (30) 

u = (v — v') (c.fEqj5|. (31) 

For our distribution, the derivatives are as follows: 


Bremsstrahlung is a significant loss mechanism for IEC 
systems; both efficiency and safety demand that our sys- 
tem does not radiate all the input energy as X-rays. 

For electron-electron bremsstrahlung we can use the 
result of Gould mm- 


3 6rem,e,e _ 20(44 37T 2 ) 3/^x2 

V 9^/tt a me y 


Xtr e n 2 e (kT e f /2 (26) 


df_ 

dv L 


and 


1 

vl 


2v l (1 + K{H)) 


f 




(32) 



where 


kef (33) 


with a the fine structure constant, T e the electron tem- 
perature. 

In the case where an ion species is in thermal equi- 
librium with the electrons, an expression for electron- ion 
bremmstrahlung power density, due once again to Gould 
is 


brem, e,i,eq 

v 
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. Ire 


llr. 


= T V ~ a ~^ neZ i niT ^ + - 24") ( 2? ) 


K(m= 2/3 exp {— /3 2 (J/ — H c ) 2 } 
y/n erfc (/? (H - H c )) 

with H = v\ + Vq + Z(f>. The term K(H) is a pure con- 
sequence of our cutoff. Without such a cutoff, it is seen 
that upscattering and heating power may be dramatically 
underestimated. 

Recalling that the average energy of a particle of 
species a is by definition 


with r = kT e /mc 2 . 

In our device, the ions and electrons are not necessar- 
ily in thermal equilibrium. Where the ions are slower 
than the electrons, the effect on bremsstrahlung is small. 
Where the ions are faster than the electrons we can sim- 
ply replace r in the above with r' = ( Tim e /T e mi)r i.e 
using the approximate electron temperature in the rest- 
frame of the ion as suggested by Jones [29: . The electron- 
ion Bremsstrahlung is computed using the larger of r and 
t'; this selection is performed at each point in space. 


(E a ) = - f f(v)E(v) d 3 v (35) 

n J 

the total heating power experienced by species a is 

-Pheat, a(z) = n a (x) ^Xf = J J dv L dvg E(v) 

(36) 

/ OO rOO 

dv e / dv L v • 3 (f a , fb,v) (37) 

-oo J 0 
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Here, use has been made of the product rule for con- 
tinuously differentiable functions: 

f du f f dv 

/ —^vd£}= / uvVidT— / u—^-d$d, (38) 

JQ OXi J v Jq OXi 

where T is the boundary of volume U. In each of our 
particular cases at least one of these terms vanishes. 

Up-scattering-driven loss power can be computed in 
the same fashion with merely a change of integration 
limits such that we integrate over that region of the dis- 
tribution holding particles that can escape and replace 
E(v) with a bound on the kinetic energy at the chamber 
wall. 


-Uipscattering loss for species a 

r°° 

— F, C ab / dv e X 
b J °° 

(vg + H c ) Vl • J(/a,/ft, v = (w a ,Vg)) 


(39) 


Note that we remain in dimensionless units; to get 
power in Watts, we multiply the above by kT (expressed 
in Joules). 


G. Power balance 


Each fusion event consumes two ions. The number 
of fusion events per unit time is Rf us = Pfus/^fus- Main- 
taining energy balance requires we replace fused ions with 
fresh ions at the pre-collision energy. Since this is neces- 
sarily less than the maximum ion kinetic energy, we can 
bound this by: 


Pion loss ^ Pion lipscattering loss + 2P { 

r 


| Qi 0coil | 

3 77' 

-C/fllS 


(40) 


t j. /o _ Useful Power Out 
^ Power in 

If the fraction of fusion power which can be captured 
is rjf then 


Useful Power Out = 77/Pfusion (41) 


Power loss in the coils will be omitted under the as- 
sumption that a real fusion power system will employ 
superconducting magnets. We will treat only the case 
where ions exert a net heating effect on the electrons; we 
will then add the electron heating power to the electron 
loss power. The continuous power input required by a 
fusion system such as has been described is at most 


Power in — -PbremT Pe lectron heating T -Pion lossT -^electron loss 

(42) 

Thus 


Q = 


dfPfusion 

Pbrem T ^electron heating T Pion loss T -^electron loss 


(43) 

The combined-cycle efficiency of a modern gas-turbine 
approaches 60% - we will thus take r] = 0.6 as a repre- 
sentative choice. 


Results 


For this paper, we restrict ourselves to low densities; in 
particular, we ensure that the particle currents are small 
enough that we can use the vacuum values for the mag- 
netic vector potential. Where this is true, we can solve 
iteratively for a self-consistent density and electrostatic 
potential structure by repeatedly recomputing densities 
and averaging the resulting 0. The electrostatic potential 
for a given density profile, coil potential, and grounded 
chamber is itself computed through an over-relaxation 
developed for this paper. When computing electron den- 
sities, the value H c for the electron population has been 
continuously adjusted such that Z(j) >= H c everywhere 
on the interior of the device and Z<p = H c at at least one 
point on the boundary. 

Electrons outside the device will be both fast-moving 
and rapidly lost to the walls. To model this, electron 
densities in the region exterior to the coils have been 
held at zero during the relaxation. 

Once the potentials Aq and (j) have been obtained ev- 
erywhere on our spatial lattice, numerically integrate is 
performed to compute the gain and loss terms introduced 


in Eqns.[37|and 39 In the process, we derive spatial maps 
of the various quantities. 

The first configuration examined is an electron-free 
system exhibiting a virtual cathode in the centre of 
the device. Parameters for this device are in table [H 
The potential-well structure with and without plasma is 
shown in Figs. 1® 0 and [8j A local virtual anode 
is formed, but the potential of this local anode is lower 
than any point on the boundary. 

The corresponding ion density is shown in Fig. [9] It 
must be understood that these 3D pictures show a slice; 
the corresponding density integrated over the azimuthal 
angle is given in Fig. [10] showing the total contribution 
from different radii. 

Azimuthally-integrated fusion density for this config- 
uration is depicted in Fisim Most of the fusion occurs 
close to the axis and near the planes of the coils. With no 
electrons to contribute to Bremsstrahlung, the dominant 
energy loss mechanism is the upscattering of ions; this is 
shown to be insignificant. The gain of the system reduces 
to 0.6 x 2.5MeV/100keV = 15 as shown in Table [TJ 

Introducing electrons allows higher densities through 
cancellation of the space-charge due to the ions. Pa- 
rameters for an example system are given in Table [TTJ 
Bremsstrahlung power is by far the dominant loss mech- 
anism in this system; one should expect the electrons in 
such a system to exhibit considerable cooling, reducing 
this loss. The proper treatment of this reduction is left 
to future work. 

The potential- well structure in F igs . [14| and [15] exhibits 
a small secondary virtual cathode - the “Poissor” struc- 
ture proposed by Hirsch [3; . The overall virtual anode is 
responsible for a great enhancement in electron confine- 
ment. Figs. [T6| and [T7| show the ion and electron densi- 
ties, respectively. The formation of a jet-like structure is 
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Parameter 

Value 

Vacuum chamber radius 

40 mm 

Vacuum chamber length 

70 mm 

Coil major radius 

20 mm 

Coil minor radius 

4 mm 

Coil separation 

35 mm 

Coil current 

700 kA-turns 

Coil potential 

-100 kV 

kT (ions) 

10 keV 

ro (ions) 

5 mm 

P (ions) 

100 

H c (ions) 

0 

Ffus 

2.5 MeV 

spatial resolution 

0.5 mm 

Computed Quantity 

Value 

Max. ion density 

4.7 x 10 16 m -3 

Total fusion rate 

2650 / s 

Total fusion power 

1.06 nW 

Ion upscattering power 

1.0 x 10 -23 W 

Q 

« 15 


TABLE I: Parameters for an electron- free system with a vir- 
tual cathode 


Vacuum Potential 



FIG. 5: Vacuum electric potential of electron- free virtual- 
cathode configuration 

notable. 

Fusion in this electron/ion configuration will occur in 
three rings, as show in Fig. [l8j For a higher-density sys- 
tem, we would alter the geometry to raise the potential 
in the ring cusp and consequently shift more of the ion 
density towards the center of the system. 

In both the cases depicted the gain, Q of the system is 
considerably greater than 1. No attempt has been made 
to optimise the parameters of the system; such a study 
awaits access to more computational power and extension 


Electric Potential 



FIG. 6: Electric potential of virtual-cathode configuration 
with ions present 



FIG. 7: Electric potential along a radial line through the 

midpoint of the device for electron- free configuration. 

of the density /potential solver to higher densities where 
simple fixed-point iteration is no longer applicable. 


Conclusions and Outlook 

It has been shown that an IEC device can achieve en- 
ergy gain if the electrodes are circular coils carrying suffi- 
cient current to prevent bombardment by and subsequent 
collection of charged particles. The required current can 
be estimated easily from the energy of the particles in 
the plasma. Gains of the order of 10 are predicted for 
a deuterium-burning system at low density. The combi- 
nation of electric and magnetic potentials in this case is 
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FIG. 8: Electric potential along an axial line through centre 

the device for electron-free configuration. 



FIG. 9: 3D slice of ion density for the electron- free virtual- 
cathode configuration 


sufficient to confine a significant electron population close 
to the device axis, admitting the possibility of a beam- 
like plasma with recirculating ions. The dimensions were 
chosen for numerical convenience - they should not be 
taken as representing limits to energy output. We see no 
reason why the device should not scale to MW outputs, 
with achievable fusion power bounded by thermal and 
mechanical properties of the materials used in construc- 
tion. 

Our results constrast with earlier work by Rider (33] , 
whose conclusions we regard as dubious, and which 
cannot be recovered unless we assume a Maxwellian, 
quasineutral plasma with the density distribution due to 
Krall (derived for non-Maxwellian distribution) and ne- 


Ion Density 



- 3 - 2-1 0 1 2 3 

Z (cm) 

FIG. 10: Ion density in integrated (r,z) space for electron- 
free virtual-cathode configuration. Black circles show the coil 
cross-section. 


Fusion Rate Densitv 



- 3 - 2-1 0 1 2 3 

Z (cm) 

FIG. 11: Fusion-rate density in integrated (r,z) space for 
electron-free virtual-cathode configuration. Black circles show 
the coil cross-section; most of the fusion is occur ing on the axis 
near the planes of the coils. 


gleet the use of the electrostatic potential for the mit- 
igation of cusp losses. We are unable to imagine an 
IEC device consistent with all of Rider’s assumptions; 
indeed concerns over self-consistency motivated much of 
our work. 

Some caveats apply; we have deliberately ignored tur- 
bulence and instabilities; we have also ignored drifts 
which can arise through inhomogeneities in the plasma 
such as the diamagnetic drift. We quote results only for 
low densities of the order of mTorr and for power out- 
puts in the nanoWatts. We have also assumed a form for 
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Parameter 

Value 

Vacuum chamber radius 

40 mm 

Vacuum chamber length 

70 mm 

Coil major radius 

16 mm 

Coil minor radius 

4 mm 

Coil separation 

28 mm 

Coil current 

500 kA-turns 

Coil potential 

-100 kV 

kT (ions, electrons) 

10 keV 

ro (ions, electrons) 

5 mm 

f3 (ions, electrons) 

100 

H c (ions) 

0 

H c (electrons) 

20.9 keV 

Ef us 

2.5 MeV 

spatial resolution 

0.25 mm 

Computed Quantity 

Value 

Max. ion density 

1.86 x 10 17 m -3 

Max. electron density 

2.04 x 10 17 m" 3 

Total fusion rate 

42360 / s 

Total fusion power 

17 nW 

Total Brem power 

0.4 nW 

Electron upscattering power 

2.7 x 10“ 4 nW 

Electron heating power 

7 x 10“ 4 nW 

Ion upscattering power 

2.7 x 10“ 9 nW 

Ion upscattering power 

1.0 x 10“ 23 W 

Q 

« 9 


TABLE II: Parameters for an electron/ion system with a 

virtual anode 


requirements to produce measurable fusion rates suggest 
the cost of building an experimental device may not be 
prohibitive, either. 
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the distribution which corresponds to a steady state of 
the Vlasov equation and thus also ignore the dynamics of 
the approach to equilibrium - the steady state of a driven 
system should not be expected to be identical to this. 

In adopting the Landau integral form for collisional 
transport, we are emphasising grazing collisions over 
strong collisions. In practice, all fusion systems must ex- 
hibit strong, inelastic, collisions - otherwise they would 
produce no fusion. The proper treatment of collisions for 
our system should involve a deeper consideration of these 
scattering events; one approach might be a monte-Carlo 
simulation involving full electrodynamics. 

Finally, we have assumed a fully- ionized system. A real 
fusion system will include neutral species; ionization and 
recombination will introduce a new loss mechanism. As 
very high levels of ionization can be achieved at the low 
densities explored in this paper, we leave the exploration 
of the neutral-species interactions for another paper. 

Future work will investigate the transition from the 
low-/? regime to the high-beta cusp-confinement regime 
treated in other papers, and develop further the theory 
of plasma confinement in charged multipole cusps. It 
would be interesting to map the parameter-space bound- 
ary wherein Q > 1. The modest volume and density 
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FIG. 12: Vacuum electric potential of electron/ion system FIG . 14: Electric potential along a radial line through the 

midpoint of the device where electrons are present. 


Electric Potential 



X < cm ) Z (cm) 

FIG. 13: Electric potential of electron/ion system with 

plasma present 



FIG. 15: Electric potential along an axial line through centre 
the device where electrons are present. 


Y (cm) Y (cm) 
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Ion Density 



Z (cm) 


FIG. 16: Ion density in the electron/ion system 


FIG. 18: Fusion rate density in the electron/ion system 


Electron Density 



— 'o 




FIG. 17: Electron density in the electron/ion system 
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Here we present the derivation of the density for our distribution away from the limit of infinitely-sharp cutoff. 


Appendix A: Density calculation 


Using the fact that 


dxxe x ! c erfc(ax 2 + b) = — -e c < e' 


b i 1 


erf (aar + b H ) + erfcfaar + b) > + C 

2ac 


dxxe x / C erfc (ax 2 + b) = ^ < erfc(6) — e^ ac + ^c^) erfc b + - — 

z [ v zac 


and therefore 


dv L F L (v L ,v e ,(f>) = ^ (erfc(/3C(vo,</>)) - exp \c{v 0 ,<j>) + -rU 1 erfc (/ 3C{v 0 ,(f> ) + 


C(v e , 4>) = v$ + Z<t> - H c = (H - H c ) 


The density integral then reduces to 


n{Ao,4>,r) = -e (ZAer 


r r oo 

|2 i e~ Z4> / dv e e -r 2 {vi+ 2 Zv e A e ) e -v 2 e erfc ^ ^2 + Z( j ) _ H 

l J —OO 

- e ^e~ H - [ dv e e - r2 (^ +2Zv ^ A ‘>) erfc (p {v 2 e + Z<p - H c ) 


We can approximate these integrals as follows: 


f(x) erfc ( g(x 2 )) dx « 2 / f(pc)dx + 


, 2g(x 2 )\ w r+ / 2g(x 2 )\ . 

1 ~^jf(x)dx+ (l y=-Jf(x)dx 


V± =M(z±) I 5(4) = 


This yields the following expression for the density: 


n(A>, </>, r) = {e" z * [J(t;_) + I(-v.) - /(-«+) + I(v+) - /(«_)] 

-e^ [J'(v'_) + I'(~v'_) ~ I'(-v' + ) + />+) - 


where 


v± ^\J±^+H a -Z</> 
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erf 


r 2 v + r 2 ZAe ^ 


— erf 


v 2 ZAq — r 2 v 


(A9) 


exp {— ( r 2 v 2 + 2 r 2 ZAov ) } 

4y / 7rf 5 

x | v^r ex P | (f 2 v + r 2 ZAoY /f 2 j erf ^ T ^ [2 f 2 (y^ — 2/3 (Z0 — i7 c )) — 2f 2 /3 — /3(2r 2 ZA<9) 2 ] 
+2r/3 [2r 2 ^ + 2r 2 ZA e \ 

(A10) 


, exp{-(rV 2 + 2r 2 ZA^')} 

1[V) ~ 4y^r 5 

x / V^rexp | (r 2 -?;' + r 2 ZA^) 2 /r 2 | erf 


r 2 v' 


+ r 2 ZAo 


) 2r 2 (V^ - 2/3 (Z</> - ff c ) - 1) - 2r 2 /3 - f3(2r 2 ZA e ) 2 ' 


+2r/3 [2rV + 2r 2 ZA 0 ] 


(All) 


J'(v f ) = — — e^ rZAe ^ [erf (rv f + vZAq) — erf ( tZAq — rv')\ 


(A12) 


and f = >/r 2 + 1. 



